UNCLASSIFIED 


AD  NUMBER 


AD878595 


NEW  LIMITATION  CHANGE 
TO 

Approved  for  public  release,  distribution 
unlimited 


FROM 

Distribution  authorized  to  U.S.  Gov't, 
agencies  and  their  contractors;  Critical 
Technology;  30  NOV  1970.  Other  requests 
shall  be  referred  to  Naval  Electronics 
Laboratory,  Center,  San  Diego,  CA  92152. 


AUTHORITY 


usnelc  notice,  26  jun  1972 


THIS  PAGE  IS  UNCLASSIFIED 


Hi 

00 

£- 

<30 

% 

BINARY  DETECTION  OF  RANDOMLY  OCCURRING  SIGNALS 
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PROBLEM 

Obtain  design  procedures  and  criteria  for  reliable  Navy  communica¬ 
tion  in  a  hostile,  an  adverse,  or  an  unstable  environment,  through  use  of 
statistical  decision  methods  and  other  signal  processing  techniques.  The  spe¬ 
cific  problem  with  which  this  report  is  concerned  is  the  detection  of  signals 
which  occur  at  random  (unknown)  times. 

RESULTS 

].  The  binary  moving-window  detector  is  a  feasible  method  of  detect¬ 
ing  signals  which  occur  at  random  times. 

2.  The  statistics  of  the  moving-window  detector  can  be  obtained  by  an 
exact  computational  method  without  resorting  to  Monte  Carlo  simulation. 

3.  The  binary  moving-window  detector  can  be  easily  implemented  for 
use  in  existing  and  future  communications  systems. 

RECOMMENDATIONS 

1.  Investigate  existing  communications  systems  and  plans  for  future 
systems  for  possible  application  of  the  binary  moving-window  detector  as  an 
inexpensive  and  effective  means  of  detecting  communications  signals. 

2.  ExterJ  the  analysis  to  cover  multilevel  quantization  of  data. 

3.  Investigate  methods  of  obtaining  analog  moving-window  detectors 
(to  avoid  theioss  due  to  quantization)  which  can  be  simply  implemented. 

4.  Analyze  the  use  of  distribution-free  statistical  procedures  in  conjunc¬ 
tion  with  moving-window  detectors. 
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INTRODUCTION 

Much  of  the  literature  on  statistical  signal  detection  is  concerned  with 
the  problem  of  choosing  between  two  hypotheses  about  a  random  variable  X 

under  the  assumption  that,  for  all  of  the  n  observations  xj, _ ,  xn,  on  X, 

exactly  one  of  the  hypotheses  is  true,  in  signal  detection  terminology,  this 
means  that  either  all  n  observations  are  of  signal-plus-noise  or  none  is.  This 
situation  presupposes  that  knowledge  is  available  as  to  when  a  signal,  if  pres¬ 
ent,  will  occur;  for  otherwise  it  would  not  be  possible  to  ensure  that  the 
above  assumption  is  satisfied.  For  example,  consider  a  communications  sys¬ 
tem  where  knowledge  is  available  at  the  receiver  as  to  when  a  message,  if 
transmitted,  will  be  received.  Then,  during  the  time  when  a  signal  can  occur, 
the  received  data  will  result  '.'mm  either  signal-plus-noise  or  noise  only  (but 
not  both).  Howeyer,  if  no  knowledge  is  available  as  to  when  a  message  will  be 
received,  the  data  will  first  result  from  noise  only  and  will  change  to  signal- 
plus-noise  at  some  unknown  time.  (This  would  be  the  ease,  for  example,  in 
signal-intercept  applications.)  It  is  this  latter  situation  which  will  be  dis¬ 
cussed  here;  that  of  detecting  signals  which  occur  at  some  unknown  times. 
Emphasis  is  placed  exclusively  on  the  detection  aspect;  that  is,  the  procedure 
is  useful  only  for  detecting  if  a  signal  Is  present  or  not  Some  other  data 
processing  is  necessary  to  determine  information  such  as  message  content. 

Discussion  will  be  confined  to  detection  methods  involving  “binary 
integration,”  a  term  used  to  describe  detection  based  on  data  which  are 
quantized  into  two  levels,  one  and  zero  (ref.  1-5).  In  particular,  a  so-called 
“moving-window  detector”  (ref.  1  and  2)  will  be  discussed,  and  methods  for 
calculating  the  detection  statistics  (false-alarm  probability,  etc.)  related  to 
this  detector  will  be  described. 

Results  will  be  given  which  will  show  how  the  detection  probability 
varies  with  SNR,  signal  duration,  and  sample  size  (the  window  “width”),  and 
a  brief  discussion  cf  the  equipment  required  for  implementation  will  be  given. 

THE  DETECTION  PROBLEM 

The  detection  problem  considered  here  is  that  of  detecting  signals 
which  appear  at  some  unknown  times  and  remain  for  a  fixed  time.  The 
problem  can  be  described  as  follows.  Let  xn,  n  =  1,  2, . . . ,  be  a  sequence  of 
independent  observations  on  the  output  X  from  a  receiver,  considered  as  a 
random  variable  with  density  function  f(x,  0),  0  an  unknown  parameter. 
Suppose  that  0  =  0O  corresponds  to  the  case  where  the  output  results  from 
noise  alone  and  that  9  =  9  j  corresponds  to  the  case  where  the  output  results 
from  signal  plus  noise  [of  some  specified  SNR,  S[0  j)j.  Suppose  also  that 
9  =  0O  for  n  =  1, 2,  —  ,  mj,  and  9  =  9 j  for  n  =  rri*  +  1, . . . ,  mj  +  L,  where 
mj  is  a  random  variable  denoting  the  observation  at  which  the  first  signal 
appears.  For  n  >  m  j  +  L,  9  =  9Q  until  the  next  appearance  of  a  signal  at  the 
random  time  m2,  etc.  It  will  be  assumed  that  two  signals  will  be  separated 
by  at  least  N  observations,  where  N  is  the  sample  size  used  in  the  test  proce¬ 
dure  (discussed  later).  It  will  also  be  assumed  that  the  duration  of  any  signal 
will  be  a  known,  fixed  number  L  of  observations. 
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As  mentioned  in  the  introduction,  only  detection  will  be  discussed 
here.  However,  a  detector  output  indicating  that  a  signal  is  present  can  be 
considered  by  the  receiver  as  an  “alarm*’  which  activates  recognition,  loca¬ 
tion  (in  time  or  frequency),  or  data  processing  stages.  Figure  1  is  a  simplified 
block  diagram  of  a  complete  system  utilizing  a  moving-window  detector  and 
a  time-lag  recorder.  The  receiver  output  is  fed  simultaneously  to  both  the 
moving-window  detector  and  the  time-lag  recorder.  If  the  moving-window 
detector  decision  is  that  a  signal  is  present,  an  alarm  causes  the  time-lag  re¬ 
corder  contents  to  be  read  out  and  processed.  The  time-lag  recorder  could 
be  a  tape  loop,  magnetic  drum,  core  storage  (for  digital  data),  or  any  conven¬ 
ient  means  of  storing  data  collected  during  a  short  time  interval.  The  actual 
length  of  time  for  which  the  time-lag  recorder  must  store  data  is  determined 
by  such  factors  as  the  expected  duration  of  signals  to  be  processed  and  the 
expected  time  required  for  the  moving-window  detector  to  produce  an  alarm 
after  the  signal  begins. 


Figure  1.  .Simplified  Block  Diagram  of  a  System  Utilizing  the  Moving-Window 
Detector  As  An  Alarm. 

THE  DETECTION  PROCEDURE 

The  method  c f  detecting  signals  (occurring  at  random  times  as  indi¬ 
cated  above)  which  will  be  discussed  here  is  the  so-called  “moving-window 
detector”  (also  called  a  moving  average  or  a  sliding  window)  (ref.  1  and  2). 
The  moving  window  of  width  N  operates  as  follows.  For  each  of  the  sub¬ 
sequences  of  observations  (xj, _ ,  xjq] ,  (x2,  •  -  - ,  *N+l)  »  —  » (xk>  -  -  *» 

xk+N-l)  >  -  -  -  >  of  length  N,  a  test  statistic  S^+N-l  (x^,  ....  *k+N-l}  *s 
computed  and,  on  the  basis  of  this  statistic,  a  decision  is  made  that  either 
0  =  0O  or  0  =  0  j.  The  test  performed  on  each  subsequence  will  be  called  a 
subtest,  and  a  specific  test  statistic  will  be  discussed  later. 

A  representative  sequence  of  observations  when  a  signal  is  present  is 


where  the  asterisk  indicates  observations  from  the  signal-plus-noise  popula¬ 
tion.  Note  that  L+N-l  consecutive  subsequences  (as  described  in  the  previous 
paragraph)  will  include  at  least  one  signal-plus-noise  observation.  This 
prompts  the  following  definition  of  a  false  alarm  and  of  a  detection:  If  no 
ooservation  in  a  particular  subsequence  is  from  the  signal-plus-noise  popula¬ 
tion  and  the  subtest  decision  is  that  0  =  0  j  (a  signal  is  present),  then  a  false 
alarm  occurs.  Likewise,  if  for  any  of  the  L+N-l  consecutive  subsequences 
for  which  at  least  one  observation  is  from  the  signal-plus-noise  population, 
the  decision  is  that  6  =  0],  then  a  detection  occurs.  Optimization  will  be 
defined  as  the  maximization  of  the  probability  of  detection,  with  the  false- 
alarm  probability  held  fixed. 

The  false-alarm  probability  will  be  denoted  by  a,  and  the  probability 
of  detection  will  be  denoted  by  D(  ),  where  the  parenthetical  entry  will 
determine  the  pertinent  variables  of  which  the  detection  probability  is  a 
function.  For  example,  in  the  context  above,  D(L,N)  denotes  the  proba¬ 
bility  of  detection  as  a  function  of  the  signal  duration  L  and  the  test  size  N. 
Formulas  for  (or  methods  for  calculating)  a  and  D(-)  will  be  given  later. 

Note  that  since  a  decision  occurs  for  each  observation,  the  average  false-alarm 
rate  (in  false  alarms  per  observation)  is  equal  to  a. 

BINARY  INTEGRATION 

A  test  which  is  of  practical  interest,  due  to  its  simplicity  and  ease  of 
implementation,  is  the  binomial  test.  When  applied  to  signal  detection,  this 
test  is  call  id  “binary  integration,”  “k  out  of  n  detection,”  “double-threshold 
detection,”  and  “coincidence  detection”  (ref.  1-6).  For  binary  integration, 
the  statistics  of  the  moving-window  detector  can  be  obtained  analytically; 
i.e.,  without  resorting  to  Monte  Carlo  simulation.  The  calculation  of  the 
statistics  of  the  binary  moving-window  detector  will  be  discussed  in  the  next 
section. 

The  binary  moving-window  detector  operates  as  follows.  The  data 
xn,  n  =  1,  2, . . . ,  are  compared  with  a  threshold  q,  and  the  random  variable 
yj  is  given  the  value  1  if  Xj  >  q,  and  is  given  the  value  0  otherwise.  Thus  the 
sequence  xn,  n  =  1,  2, . . . ,  is  converted  into  the  binary  sequence  yn,  n  =  1 , 
2,...  .  At  the  k**1  observation  the  statistic 

k 

Sk-  I  y\  <0 

i=k-N+l 

is  calculated  and  compared  with  a  threshold  V  and,  if  Sj<  >  V,  the  decision 
is  made  that  0  =  0i  (a  signal  is  present);  otherwise,  a  decision  is  made  that 
0  =  0O* 


*Some  authors  define  a  detection  as  the  occurrence  of  a  threshold  crossing;  that  is,  a  sig¬ 
nal  decision  occurs  at  the  k^1  observation  only  if  Sj^  =  V  jnd  Sjc_j<  V.  (Ref.  7,  page  767.) 


Under  the  previous  assumptions  about  the  density  of  X.  it  follows 
that 


p(yn=n o  =  o0)  =  J  f(x,o0)dx 
q 

CO 

p  (yn=ll0  =  0l)  =  /.f(x,Cjjdx 

q 


(2) 


(3) 


Denoting  these  probabilities  by  p0  and  pj,  respectively*  it  follows  that* 


(4) 


The  following  procedure  is  used  to  optimize  the  test,  in  the  sense  of  maxi¬ 
mizing  D(V,N)  forgiven  L.  a,  and  SNR  (ref.  6).  (Here,  D(V,N)  denotes  the 
probability  of  detection  as  a  function  of  the  threshold  V  and  the  test  size  N.) 
Let  N  be  fixed.  For  each  V,  V  =  1,  2, . . .  ,FI,  equation  4  is  solved  for  the 
value  p0(V)  for  which  the  given  false-alarm  probability  a  is  obtained".  For 
each  p0(V),  equation  2  is  solved  for  the  threshold  q(V),  and  for  each  q(V), 
equaOion  3  is  solved  for  pj(V).  The  pairs  (po(V),  P](V)),  V  =  1 ,  2, . . . ,  N, 
so  obtained  are  used  (in  a  method  to  be  described  in  the  next  section)  to 
calculate  D(V,N),  and  the  value  of  V  for  which  D(V,N)  is  maximum  is 
selected.  Now  N  is  allowed  to  vary  and  the  above  procedure  is  followed  for 
each  N.**  The  test  size  N  is  then  chosen  as  the  value  of  N  for  which  D(V,N) 
is  maxirtiized,  having  chosen  the  optimum  value  of  V  for  each  N.  Results 
will  be  given  later  to  show  how  the  probability  of  detection  varies  with  N, 
for  fixed  L  and  a. 


*Equation4is  valid  for  any  particular  subtest,  and  the  average  false-alarm  rate  is  a  false 
alarms  per  observation  even  though  the  subtests  are  not  independent. 

**It  will  be  seen  from  results  presented  later  that  usually  only  values  of  N  in  the  neighbor¬ 
hood  of  L  need  be  tried  in  the  optimization  procedure;  that  is,  N0pt  ~  L.  Also,  for 
many  of  the  probability  densities  of  interest,  the  optimum  threshold  V  (for  fixed  N) 
can  be  determined  without  need  of  considering  all  possible  values  for  V.  This  is  the 
case,  for  example,  when  the  detection  probability  D(V)  (for  fixed  L,  N,  and  a)  is  uni- 
modal.  It  is  then  only  necessary  to  determine  the  first  value  of  V  for  which  D(V  + 1) 

<  D(V).  Finally,  it  is  sometimes  possible  to  solve  equations  2  and  3  for  pj  as  a  function 
of  p0,  and  thus  avoid  the  intermediate  step  of  determining  q  (ref.  6,  page  349). 


CALCULATING  THE  PROBABILITY  OF  DETECTION 

The  method  used  for  calculating  the  probability  of  detection  for  the 
binary  moving-window  detector  is  an  application  of  the  theory  of  finite 
Markov  processes*  (ref.  8-10).  This  theory  cannot  be  applied  directly  to  Sj- 
since  {Sj^}  ><=]  is  not  a  Markov  process  (ref.  8).  However,  if  all  possible 
configurations  of  ones  and  zeros  in  the  k1*1  subsequence  are  considered,  then 
finite  Markov  process  theory  can  be  applied;  as  discussed  below. 

The  binary  moving-window  can  be  viewed  as  follows.  Suppose  that 
after  each  observation,  the  contents  of  an  N-bit  shift  register  are  shifted  left 
by  one  bit  and  the  binary  datum  for  that  observation  is  inserted  in  the  Tight- 
most  bit.  Then,  after  the  k^1  observation,  the  content  of  the  shift  register 
is  the  k^1  binary  subsequence,  and  Sjj  is  the  sum  of  the  I's  in  the  shift  regis¬ 
ter.  If  we  let  be  the  number  obtained  by  regarding  the  shift  register  con¬ 
tent  as  a  binary  number,  then  {Rj.|  is  a  finite  Markov  process  with 
possible  states  0, 1,  —  ,2^-1.  That  jRj.}  is  a  Markov  process  follows  from 
the  fact  that  the  possible  values  of  Rjj  are  determined  completely  by  the 
value  of  R^-i  and  yj..  Note  that  a  left  shift  of  an  N-bit  shift  register  is  equiv¬ 
alent  to  multiplication  oy  two.  modulo  2^_  Hence,  if  Rj-_j  -  i,  then  Rj.  = 
(2i)mod  2^  if  yj.  =  0  and  R^  =  1  +  (2i)mod  2^,  if  y^  =  1.  1 '  p(n)  is  the 
probability  of  a  1  at  the  ri^1  observation,  then  the  transition  probabilities  are 


jiRn-1 


1  -p(n)  if  j  =  (2i)mod  2^ 
p(n)  if  j  =  1  +  (2i)mod 
0  otherwise 


(5) 


The  transition  matrix,  for  the  states  ordered  as  the  numbers  0,  1 , ....  2^  -1 . 
has  the  form 


q(n) 

0 


p(n) 

0 

/ 


0 

q(n). 


0 

p(n) 


0  0 
0  0 


P(n)  = 


0  0  0  0 

q(n)  p(n)  0  0 

0  0  q(n)  p(n) 


q(n) 

0 

0 


P(n) 

0 

0 


(6) 


0  0  0 


q(n)  p(n) 


where  q(n)  =  1  -  p(n).  Notice  that  the  lower  half  of  the  matrix  is  identical 
with  the  upper  half  since  Rn  is  independent  of  the  high-order  bit  of  Rn_] . 


*A  finite  Markov  process  Is  a  stochastic  process,  which  can  take  on  only  a  finite  numbei 
of  values,  such  that  the  outcome  at  the  k^1  stage  of  the  process  depends  on  the  outcome 
at  the  (k-l)st  stage  and  on  r.o  previous  outcome.  If,  additionally,  the  outcome  at  the 
k*k  stage  is  independent  of  k,  the  process  is  called  a  Markov  chain. 


To  avoid  cumbersome  subscripts,  in  the  following  we  will  let  W  =  2^. 

if  the  components  of  the  probability  vector  IJn=  .^jtqC11), _ ,  T\y_](r)) 

denote  the  probability  that  Rn  lias  values  0, 1, . . . ,  "W-l,  then  the  corre¬ 
sponding  probability  vector  Iin+-j  can  be  found  from 

^n+1  =  nnP(n+l )  (7) 

F rom  the  form  of  P^n)  in  equation  6,  note  that  the  indicated  matrix  multipli¬ 
cation  in  equation  7  reduces  to 

^  =  +  Trj4\y. /2^3  5(n+^ ) 

troj+t  (n+^  )  —  [jr^11)  -r  7rj4--\yy2^n^3  P(n^ )  (8) 


For  the  problem  treated  here,  it  is  assumed  that  the  probability  of 
entering  a  1  into  the  window  in  the  norsignal  case  has  a  constant  value  p0, 
independent  of  n.  It  will  be  shown  that  after  N  observations  tlie  process 

reaches  equilibrium  in  the  sense  that  31^^ = IIjvj,  n  =  1,  2, _ ,  until  the 

process  is  disturbed  by  the  occurrence  of  a  signal.  Specifically,  ifKj  is  de¬ 
fined  as  the  number  of  T’sin  the  binary  representation  of  the  number  j,  then 

I^N+n "  ^0’ '  ■  * 1  ^w-l)  ■ 

®  P0Kj  (3  "  Po)  N_Kj;  j  =  0, 1 , . . -  ,  W-3  (5b) 

To  see  that  equation  9  is  trueforn=  0,  note  that  after  N  observations  the 
content  of  the  shift  register  represents  the  results  ofN  independent 
Bernoulli  trials  the  ordered  Tesults  of  which  have  probabilities  given  by 
equation  9b.  Next,  note  that  Kj+-\y/2~  Kj+1  -  K?j+j  fori^  0, 1, ... , 

W 

~-l.  Using  this  and  equation  9b  in  equation  ^,  we  have 

*2i(N+,)  -  [poKi(l-Po)N',Ki+I)]  (1-Po) 

=  P.Ki(l-Po)Mi-"2P 

'2i+l*N+1)  =  PoK>(l-Po)N-|K-+"  Po 

=  PoKl+‘  (l-P0)N'fKi+,:>  'W* 

which  implies  that  equation  9  is  true  for  n=  1,  2, ...  , 

As  stated  previously,  a  signal-present  decision  occurs  whenever  the 
number  of  ones  in  the  window  is  greater  than  or  equal  to  the  threshold  V- 


Equivalently,  a  signal-present  decision  occurs  whenever  the  outcome  Rk  is  a 
member  of  the  subset  D  of  states  such  that  jeD  if  Kj  >  V.*  (Kj  is  the  num¬ 
ber  of  1  ’s  in  the  binary  representation  of  the  number  j.)  A  state  in  D  will  be 
called  a  decision  state. 

From  the  definition  of  a  detection,  when  a  signal  is  present  there  are 
exactly  N+L-l  consecutive  observations  for  which  a  correct  signal-present 
decision  (not  a  false  alarm)  can  occur.  The  probability  of  detection  defined 
previously  can  be  restated  as 

D(L,N)  =  1  -P(Rk^  D  for  k  =  m+ 1 . m+L+N- 1)  (10) 

where  it  is  assumed  that  the  first  signal  observation  is  xm+j .  In  words,  equa¬ 
tion  10  states  that  1-D(L,N)  is  the  probability  that  the  process  fails  to  enter 
a  decision  state  as  the  signal  data  “passes  through”  the  window.  Calculating 
the  probability  of  detection  is  accomplished  by  the  following  recursive  pro¬ 
cedure. 

Assume  that  there  have  been  at  least  N  observations  of  noise  only,  so 
that  the  probability  vector  H0  has  the  form  given  by  equation  9.  (For  con¬ 
venience  of  notation,  the  observations  have  been  reindexed  so  that  the  first 
signal  observation  is  x j.)  Let  Aq=  (  a.q^),  . . .  where Aj/01  = 

7Tk^)  if  k^D,  and  Xk(^)  =  0  otherwise.  Compute  A  j  =  (AqI  'I.  . . .  .  A\y_j( '  l) 

from  A0  by  'first  using  equation  8  and  then  setting  Aj,^)  =  0  if  keD.  (That 
is.  A]  =  AqPO)  with  all  components  of  A-j  corresponding  to  decision  states 
set  equal  to  zero.)  Continue  computing  An+j  from  An  in  this  fashion  until 
AN+L-1  *s  °bta‘ne(i-  In  these  computations,  p(n)  (the  probability  of  entering 
a  1  into  the  window)  is  given  by  p(n)=  p  j  for  n  =  1 , 2, .  . . ,  L,  and  by  p(n)  = 
p0  for  n  =  L+l , . . . , 'N+L-l ;  corresponding  to  the  situation  where  the  signal 
data  is '‘entering”  the  window  and  “leaving”  the  window,  respectively.  The 
probability  of  detection  is  given  by 

D(L,"N)  =  1  -  y  Aj(N+L- 1 )  ( l  1 ) 

jfD 

That  is.  l-iDfL.'N)  is  the  sum  of  the  components  of  1  corresponding 

to  nondecision  states 

The  validity  of  equation  11  will  be  proved  by  showing  that,  for  n  =  1 . 

2. . . . , N+L-l. IP [!Rri£D.iRn_!#D, . . .  ,R0«®]  =  J  Xj(n),  with  Aj(»)  obtained 

j  4D 

as  indicated  above.  Asbefore,'W=  2^.  Note  that 


*iThemotation  jeD  means  “j  isa  member  of  the  set  D.”  Similarly,  te'D  means  “j  is  not  a 
member  of  the  sct!D.” 


(12) 


F  [Rj  =j,  Rq^D]  =£  P[R,  =jJR0  =  i]P[R0  =  i] 
i<?D 

W-! 

=  ^  P[Rl  =j|R0  =  i]Xj(0)  H3) 

H> 

W-I 

=  £  Plj(l)Xi<0>  (14) 

i=0 

and  (by  definition)  Xj(  1 )  =  P[R]  =j,  Rq^D]  if  j^D  and  XjO)  =  0  if  jeD. 
Equation  13  follows  from  equation  12  since  Xj(0)  =  P  [Rg  =  ij  if  iVD  and 
Xj(0)  =  0  if  ieD,  and  equation  14  follows  from  equation  13  by  the  definition 
of  pjj(l)  in  (5).  Hence,  P[Rj^D,  Rq^D]  =  x/^.  Now,  suppose  x/m^  ~ 

#D 

P[Rm  =  i,  Rm_i^D, . . . ,  Rq/D]  if  \4D  and  Xj(m)  =  0  if  ieD,  i  =  0 . W-l. 

Then 


p  [Rm+1  "J*  Rm^D>  ’  R0^D] 

=  J  P[Rm+l=jlRm=  i,Rm^r.,...,Ro0] 
tfD 

•P[Rm  =  i,Rm_^D,-..,Ro^]  CI5) 

Applying  the  Markov  property  (see  the  footnote  on  page  7)  to  the  conditional 
probability  in  equation  15  and  using  the  definitions  of  Xj(m)  and  pjj(m+l), 
we  have 


W-l 

p[Rm+l  =J’  Rm^D>  •  •  • .  Ro^D]  ~  ^  p[Rm+l  =jlRm 

i=0 
W-l 

=  Pjj(m+J  )Xj(m; 
i=0 

the  right-hand  side  of  which  is  an  application  of  equation  8  to  Am.  Setting 
Xj(,n+1)  =  P[Rm+i  =-j,  Rm^D, . . . ,  Rq^D]  if  HD  and  Xj(m+1)  =  0  if  jeD,  it 

follows  that  P  [Rm+]^D, . . . ,  Rq^P]  =  ^  X^m+1^  and  equation  1 1  follows 

i<fD 


=  i]Xi(m) 
(16) 

07) 


by  induction. 


RESULTS 

Results  are  given  for  three  probability  densities,  the  Rice,  Rayleigh, 
and  Gaussian  (normal)  densities.  The  Rice  density  is  used  to  describe  the 
output  X  from  a  linear  envelope  detector  when  the  SNR  is  constant  (non¬ 
fluctuating).  This  density  is  given  by 

f^x,0^  =  x  expj-^p  +  0^j  I0 (xn/201) dx,  x^O  (18) 

where  0  j  is  the  SNR,  and  I0(-)  is  the  modified  Bessel  function  of  the  first 
kind,  order  zero  (ref.  i  1).  The  corresponding  no-signal  density  is  given  by 

f(x, 0)  =  x  exp[-x-/2] ,  x>0  (19) 

which  is  obtained  from  equation  1 8  by  setting  0  j  =  0  =  0O.  Solving  (2)  for 
q,  given  pG,  results  in 

q=(-2  1nPo  f  (20) 

The  solution  to  equation  3  for  pj,  given  pQ,  is  obtained  from  the  so-called 
Q-function* 

Pj  =  Q  (y/2dj,  q)  =  /  x  exp ^  +  B  l)]  Io(x\/20T)dx  (21) 

q 

which  is  tabulated  (ref.  12),  or  can  be  calculated  recursively  (ref.  13). 

The  Rayleigh  probability  density  can  be  used  to  describe  the  output 
Y  from  a  linear  envelope  detector  when  the  SNR  fluctuates  according  to  the 
Swerling  case-2  (ref.  14).  This  density  is  given  by 

f(y>8i)=i7rlexp 

where,  in  this  case,  0  j  is  the  average  SNR,  averaged  over  all  possible  fluctua¬ 
tions.  The  corresponding  no-signal  distribution  is  again  given  by  equation  19. 
The  probability  pj  can  be  obtained  in  terms  of  p0  directly,  without  the 
intermediate  step  of  calculating  q.  Since  the  solution  to  equation  2,  using 
_n2 12 

equation  19,  is  p0  =  e^  and  the  solution  to  equation  3,  using  equation 
22,  is  pj  =  e-^2  ( 1+0  0,  it  follows  that 

P,  =  Po1/(,+(’1)  <23) 


L  2(1+"i) 


y  >0 


(22) 


oo 

*The  Q-function  is  defined  by  Q(a,  t)  =  j  x  exp  [-(x2  +  a2)/'/]  I0(ax)  dx. 

t 
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It  is  noted  that  the  results  to  be  given  for  the  random  variable  X  (and 
Y)  described  above  are  valid  also  for  any  random  variable  W  =  T(X)  (Z  =  T(Y)), 
when  T  is  a  one-to-one  monotone  increasing  transformation.  Hence,  the 
results  obtained  utilizing  densities  associated  with  a  linear  detector  are  also 
valid  for  a  square-law  detector,  with  the  transformation  T  defined  by  T(u)  = 
U"/2. 

The  Gaussian  density  can  describe  the  receiver  output  under  a  variety 
of  situations.  When  a  signal  is  present,  it  is  assumed  that  the  associated  proba¬ 
bility*  density  is 

f(x,  0  j)  =  (2tt)-,/2  exp  f-  (x-0  ij  2/2j  (24) 

where  the  SNR  is  0 \~I2.  When  no  signal  is  present,  the  density  is 

f(x,  0)  =  (2jr)-^exp  [-x2/2]  (25) 


The  equation  relating  p  j  and  pQ  (used  in  lieu  of  equations  2  and  3)  is 


Pj  =  erfc 


where 


(26) 


OO 

crfc(q)  =  (2ff)->/' /c-t2/2dt  (27) 

q 

Figure  2  shows  the  probability  of  detection  D(V)  as  a  function  of 
the  threshold  V  for  several  values  of  the  SNR,  with  a,  N,  and  L  fixed.  For 
each  curve,  the  circled  point  indicates  the  maximum  value  of  D(V),  and  the 
corresponding  optimum  value  of  V  can  be  read  from  the  abscissa.  These 
results  indicate  that  the  optimum  value  for  V  is  relatively  insensitive  to  SNR. 
A  similar  result  was  presented  in  Worley  (ref.  6).  Although  insensitive  to  SNR 
(for  each  probability  density),  the  results  of  figure  2  indicate  that  knowledge 
of  the  probability  density  is  necessary.  For  example,  using  V  =  4,  which  is 
optimum  for  the  Rician  model,  significantly  degrades  D(V)  if  used  for  the 
Rayleigh  model.  It  is  expected  that  for  large  values  of  N  (and  L),  the  dif¬ 
ference  in  the  optimum  values  of  V  for  the  three  models  will  become  more 
pronounced.  That  is,  as  N  increases  it  becomes  increasingly  more  important 
that  the  probability  density  describing  the  receiver  output  be  known,  so 
that  the  optimum  value  of  V  can  be  obtained.  No  results  are  available  at 
present  to  support  this  claim;  however,  related  results  are  presented  in  Worley 
(ref.  6). 
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Figure  2.  Probability  of  detection  D(V)  as  a  function  of  the  threshold  V,  for 
several  values  of  the  SNR. 

Figure  3  shows  the  probability  of  detection  D(L,N)  as  a  function  of 
the  test  size  N,  for  several  values  of  the  signal  duration  L  where,  for  each  N 
and  L,  the  optimum  V  is  chosen.  Notice  that  for  both  the  Gaussian  and 
Rician  models,  D(L,N)  has  a  distinct  peak  which  occurs  for  L=N.  However, 
for  the  Rayleigh  model,  D(L,N)  is  relatively  flat  (as  a  function  of  N  for  each 
L),  with  a  maximum  at  N=l.  These  results  again  indicate  the  importance  of 
knowing  the  probability  density  associated  with  the  receiver  output  if 
optimization  is  to  be  achieved. 
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IMPLEMENTATION  OF  THE  BINARY  MOVING-WINDOW  DETECTOR 

As  mentioned  earlier,  binary  integration  is  of  interest  due  to  the  ease 
in  which  a  system  for  binary  integration  can  be  implemented.  One  very  im¬ 
portant  application  of  the  binary  moving-window  detector  is  the  detection  of 
communication  signals  (from  both  friendly  and  unfriendly  sources)  which 
occur  at  random  times.  If  the  number  of  data  samples  which  must  be  col¬ 
lected  and  stored  is  large,  then  binary  integration  may  well  be  the  only  feas¬ 
ible  method  of  detection  due  to  the  speed  at  which  computations  must  be 
made.  A  brief  description  of  a  particular  method  of  implementing  the  binary 
moving-window  detector  for  detecting  communication  signals  will  be  given. 

A  similar  and  more  complete  description  (in  terms  of  a  multiple-range-bin 
radar  system)  is  found  in  Dillard  (ref.  2). 

A  block  diagram  of  the  binary  moving-window  detector  of  width  L. 
as  applied  to  the  detection  of  randomly  occurring  communication  signals,  is 
shown  in  figure  4.  The  inputs  labelled  S,  S,  and  C  are  illustrated  in  figure  5. 
The  input  S  (SAMPLE  PULSES)  is  a  sequence  of  pulses  occurring  at  equally 
spaced  intervals  corresponding  to  the  rate  at  which  the  receiver  output  is 
sampled  and  quantized.  The  input  S  is  the  logical  complement  of  S  when 
considered  as  a  logical  input  to  and-gate  2.  (In  figure  5,  relatively  negative 
voltages  are  considered  logical  “ones”  and  relatively  positive  voltages  are 
considered  logical  “zeros.”)  The  input  C  (CLOCK  PULSES)  is  a  sequence  of 
pulses  occurring  at  L  times  the  rate  of  the  SAMPLE  PULSES,  where  L  is  the 
window  width. 

The  heart  of  the  detector  is  the  delay  line,  which  is  a  clocked  (digital) 
recirculating  memory.  The  total  delay  is  adjusted  to  exactly  L+l  CLOCK 
PULSE  periods,  where  L  is  the  window  width.  Logically,  the  delay  line  can 
be  considered  as  an  (L+l)-stage  shift  register  which  shifts  at  the  clock  rate.* 

At  each  SAMPLE  PULSE  the  receiver  output  X  is  compared  with  the 
reference  q  and  if  X  >  q  the  quantizer  output  Q  =  1 ;  otherwise  Q  =  0.  _Also, 
at  the  occurrence  of  each  SAMPLE  PULSE,  and-gate  2  is  disabled  by  S,  and 
thus  inhibits  recirculation  of  the  corresponding  data  bit.  And-gate  1 ,  being 
enabled  by  the  SAMPLE  PULSE,  causes  a  1  or  0  to  be  “written”  into  the 
delay  line  accordingly  as  Q  =  1  or  Q  =  0.  Thus,  the  old  data  which  were  pre¬ 
vented  from  recirculating  are  replaced  with  new  data  (from  Q,  via  and-gate  1). 
Between  occurrences  of  SAMPLE  PULSES,  and-gate  1  is  disabled  (S  is  at 
logical  zero)  ar.d  and-gate  2  is  enabled,  allowing  the  data  to  recirculate.  Since 
the  delay  line  has  total  delay  of  exactly  L+l  CLOCK  PULSE  periods,  a  data 
bit  which  is  inserted  after  the  j^1  SAMPLE  PULSE  will  precess  by  one  CLOCK 
PULSE  period  for  each  trip  through  the  delay  line.  After  a  total  of  L  trips, 
the  data  bit  will  have  precessed  to  the  point  where  it  appears  at  and-gate  2 
when  S  =  0  and  is  thus  replaced  by  new  data.  Therefore,  the  input  to  the 
delay  line  between  the  occurrence  of  SAMPLE  PULSES  is  the  binary  data 
from  each  of  the  previous  L  receiver  quantizations. 


*ln  many  applications  it  is  more  convenient  and  economical  to  actually  use  a  shift  register 
rather  than  a  (digital)  delay  line.  However,  when  L  is  large,  the  use  of  a  shift  register  may 
be  much  more  costly  than,  for  example,  the  use  of  a  magnetostrictive  delay  line. 
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Figure  4.  Block  Diagram  of  the  Moving-Window  Detector. 


SAMPLE  PULSES  (S)  -[J - V 


CLOCK  PULSES  (Ct-y-ypy-^j^  - - - JIJTJTJTT' 

K - L  PULSES - - - ^ 


DELAY  LINE  K - TOTAL  DELAY  -  L  + 1  CLOCK  PULSES - >j 

Time-^- 

Figure  5.  Timing  Relationships  Among  the  Sample  Pulses, 

Clock  Pulses,  and  the  Total  Delay  of  the  Delay  Line. 

Counter  S  is  reset  by  each  SAMPLE  PULSE.  The  input  to  the  delay 
line  together  with  the  CLOCK  PULSES  are  applied  to  and-gate  3.  As  a 
result,  counter  S  counts  the  total  number  of  “ones”  obtained  from  the  L 
previous  quantizations  of  the  receiver  output.  That  is,  counter  S  computes 

k 

Sk=  2  * 

i=k-L+l 

where  yj  is  the  binary  datum  from  the  i^1  quantization  of  the  receiver  output. 
The  contents  of  counter  S  are  compared  with  the  threshold  V  and,  if  S^  >  V, 
a  signal-present  decision  occurs.  Thus,  after  each  receiver-output  quantiza¬ 
tion,  the  comparator  indicates  if  the  threshold  V  is  exceeded. 
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CONCLUSIONS 

The  detection  of  randomly  occurring  signals  can  be  accomplished  by 
a  binary  moving-window  detector.  The  detector  is  easily  implemented  and 
can  be  applied  to  existing  and  future  communications  systems,  especially  for 
use  as  an  “alarm”  to  initiate  location,  recognition,  or  data  processing  stages. 

The  statistics  of  the  binary  moving-window  detector  can  be  obtained 
by  a  direct  calculation,  without  resorting  to  Monte  Carlo  simulation.  Methods 
for  obtaining  these  statistics  are  easily  programmed  for  computer  calculation. 

RECOMMENDATIONS 

1.  Investigate  existing  communications  systems  and  plans  for  future 
systems  for  possible  application  of  the  binary  moving-window  detector  as  an 
inexpensive  and  effective  means  of  detecting  communications  signals. 

2.  Extend  the  analysis  to  cover  multilevel  quantization  of  data. 

3.  Investigate  methods  of  obtaining  analog  moving-window  detectors 
(to  avoid  the  loss  due  to  quantization)  which  can  be  simply  implemented. 

4.  Analyze  the  use  of  distribution-free  statistical  procedures  in  con¬ 
junction  with  moving-window  detectors. 
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IS,  ABSTRACT 


Use  of  binary  moving-window  techniques  for  detection  of  randomly  wCu7iiii^  SiCmuiS  '<nuS  ir.vesti- 
gated.  Statistical  decision  methods  appear  to  offer  a  means  for  implementing  a  reliable  signal 
detection  and  processing  system  in  a  hostile,  adverse,  or  unstable  environment.  The  statistics  of 
the  moving-window  detector  can  be  obtained  by  direct  computation,  without  resorting  to  ‘Monte 
Carlo’  simulation;  and  are  easily  programmed  for  computer  calculation. 
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